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ABSTRACT 

The intermediate Helium subdwarf B star LS IV-14° 116 is a unique object showing extremely peculiar atmospheric abundances as 
well as long-period pulsations that cannot be explained in terms of the usual opacity mechanism. One hypothesis invoked was that a 
strong magnetic field may be responsible. We discredit this possibility on the basis of FORS2 spectro-polarimetry, which allows us to 
rule out a mean longitudinal magnetic field down to 300 G. 

Using the same data, we derive the atmospheric parameters for LS IV-14°116 to be T eff = 35,150±111 K, logg = 5.88±0.02 and 
log lV(He)/N(H) = -0.62±0.01. The high surface gravity in particular is at odds with the theory that LS IV-14°116 has not yet settled 
onto the Helium Main Sequence, and that the pulsations are excited by an e mechanism acting on the Helium-burning shells present 
after the main Helium flash. 

Archival UVES spectroscopy reveals LS IV-14°116 to have a radial velocity of 149.1 ±2.1 km/s. Running a full kinematic analysis, 
we find that it is on a retrograde orbit around the Galactic centre, with a Galactic radial velocity component t/=13.23±8.28 km/s and 
a Galactic rotational velocity component V=-55.56±22.13 km/s. This implies that LS IV-14°116 belongs to the halo population, an 
intriguing discovery. 
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atmospheres — Stars: oscillations — Stars: chemically peculiar 


1. Introduction 

Pulsating subdwarf B (sdB) stars have received a lot of attention 
over the last years. Having become one of the success stories 
of asteroseismology, they hold the key to a more mature under¬ 
standing of the formation and evolution of hot subdwarfs. While 
it is commonly accepted that sdB stars are compact, evolved, 
He-core burning objects that lost too much of their H-envelope 
before or at the He-flash to sustain H-shell burning, the details 
of their evolution and in particular the mass loss mechanism re¬ 
main unclear. A large fraction of field subdwarf B stars are found 
to reside in binary systems and may have been stripped of their 
envelope mass by binary interactions involving Roche lobe over¬ 
flow and/or a common envelope phase dHan et al.l l2002[j2003|h 
poten tially resulting in a late Helium-flash fe.g.. iBrown et alJ 
1200 ill . Single sdB stars have been explained by mergers involv- 
ing at least one white dwarf (lHanll2008t IClausen & WadelliZOl ill 
or a Helium-enriched subpopulation in Globular Clusters (e.g., 
iD’Antona et alJl2005h . 

First discovered in the late nineties dKilkennv et al.lfl997tl . 
the rapidly pulsating sdB stars show luminosity variations on 
timescales of 100-200 s and are found in a well-defined in¬ 
stability strip between ~29,000-36,000 K. The fast pulsations 
have been modelled very successfully in terms of pressure ( p )- 
mode instabilities excited by an opacity (/©mechanism associ- 

* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (proposal ID 
093.D-0680 and 087.D-0950) 


ated with a local overabund ance of iron in the driving region 
dCharoinet et alJH~996l.ll 997 1. While sdB stars are overall metal- 
poor, diffusion (the competitive action between radiative levi¬ 
tation and gravitational settling) creates a non-uniform abun¬ 
dance profile of heavy elements as a function of depth. As a 
consequence, all sdB stars observed so far show chemically pe¬ 
culiar atmospheric abundances, and the vast majority of them 
(—90%) are Helium-deficient (ranging from no detected Helium 
to nearly Solar). Diffusion is also the key ingredient to under¬ 
standing the pulsation driving in a well-defined instability strip 
in surface gravity - effective temperature (log g - 7/n j space: de¬ 
pending on the exact balance between the levitation of heavy el¬ 
ements caused by the radiation pressure (dependent on T e g ) and 
the gravitational settling (dependent on logg) an iron-”bump” 
is formed, causing a sharp peak in opacity at the right depth to 
drive pulsations, or not. The depth at which pulsations are ex¬ 
cited depends critically on their type: the rapid p-mode pulsa¬ 
tions are driven in the outer envelope regions, whereas in con¬ 
trast the longer-period gravity modes (g-modes) are produced 
deeper inside the star. And indeed, the iron opacity mechanism 
also gives rise to a second type of sdB pulsator with longer 
(P ~2000-8000 s) g-mode pulsations, members of which are 
located in a cooler instability strip between 22,000-29,000 K 
dGreen et alj|2003i) . It is then quite fitting that stars located at 
the intersection between the p- and g-mode instability regions 
appear to be hybrids and show both rapid and slow luminos¬ 
ity variations (iSchuh et al.li2006l ). The two types of sdB pulsator 
have now been very successfully analyzed using asteroseismol- 
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ogy (e.g. JVan Grootel et alll2013i : ICharpinet et al.ll201 I I. where 
the internal parameters of the star are derived from the observed 
pulsation spectrum to a high accuracy and can be used to start 
constraining the proposed evolutionary scenarios. 

Given that the pulsation properties of the two types of sdB 
pulsator are so well defined, it was a surprise when long-period 
(P - 2000-5000 s), multi-periodic luminosity variations were 
detected in the hot (T P « -35,000 K) sdB star LS IV-14°116 
(lAhmad & Jeffervl 12005c iGreen etafl 1201 lh . While exhibiting 
pulsation periods typical of the g-mode pulsators, this unique 
star is located right in the middle of the hotter p-mode instabil¬ 
ity strip (see Fig. 5 of [Green et alJl201~Th . where the longer pul¬ 
sations observed cannot be explained by the ^-mechanism. As 
an alternative, it was suggested that the observed pulsations may 
be excited by an e-mechanism acting in the He-burning shells 
tha t app ear before the star settles in the core He-burning stage 
(iMiller Bertolami et al.ll2()TTh . 

Compared to the other, H-rich sdB pulsators, LS IV-14°116 
is also different in that it shows a mild atmospheric Helium- 
enhancement (log A(He)/N(H) - -0.6). This makes it an inter¬ 
mediate He-rich sdB, and as such a potential asteroseismic probe 
of a different evolutionary phase to the other sdB pulsators. For 
instance, it has been suggested that He-sdBs are the immediate 
progenitors of H-sdBs, before the Helium has had time to be de¬ 
pleted from the atmosphere due to gravitational settling ( Jeffery 
[2012 ). Apart from being the only He-rich sdB pulsator, LS IV- 
14° 116 also shows a very strange abundance pattern even for 
a subdwarf B star, with Ge, Sr, Y and Zr enhanced by factors 
of up to 10,000 compared to Solar (iNaslim et al.112011 ). One 
of the most plausible explanations invoked for these observed 
overabundances was radiative levitation in a particularly quiet 
atmosphere invoked by a strong (> kG) magnetic field. Such a 
magnetic field would presumably affect the driving of pulsation 
modes, and potentially shift the g-mode instability strip to higher 
temperatures, thus naturally explaining the pulsations observed 
in LS IV-14°116 . 

The main aim of the study presented here was to check LS 
IV-14°116 for the presence of a magnetic field with a disk- 
averaged longitudinal component of a few hundred Gauss or 
higher using the spectro-polarimetric capability of FORS2 at 
the VLT on Cerro Paranal, Chile. This facility has over the last 
decade been used extensively to study stellar magnetism, and 
led to claims of magnetic field detections in a variety of stars 
across the Hertzsprung-Russell diagram, including hot subd¬ 
warfs (iO’Toole et al.l2005l ). However, several of these detections 
(including those presented for sdB stars) were later shown to be 
artefacts of the data reduction procedure in a systematic analy¬ 
sis of FORS2 spectro-polarim etric archival data (Bagnulo et al. 
120121: lLandstreet et al.ll2012ah . An independent analysis based 
on spectro -polari metry obtained with ESPaDOnS at the CFHT 
([Petit et a~L!l2Q 12 ) also failed to measure a magnetic field in two 
hot subdwarfs for which a positive detection had previously been 
claimed. In this context, it was with great interest that we noted 
the first detection of a very strong (several hundred kG) magnetic 
field in a He-rich subdwarf O star dHeber et al . 2013?) . Given that 
magnetic fields are relatively common among white dwarfs, they 
should also be present for a sizeable fraction of their hot subd¬ 
warf (sdB and sdO) progenitors. With its puzzling pulsational 
properties and atmospheric abundance profile, LS IV-14°116 
was then a prime candidate for the first magnetic sdB star. 

In the following sections we present the FORS2 spectro- 
polarimetry gathered and the resulting null detection of a mag¬ 
netic field. We then move on to the exploitation of the FORS2 
total intensity spectra and supplement these with UVES archival 


data to show that LS IV-14°116 is very likely a halo star. The pa¬ 
per ends with a discussion of the open questions regarding this 
unique object. 


2. No evidence for a magnetic field 

Our spectro-polarimetric data were gathered with FORS2/VLT 
in the second half of the night on 4th and 5th June 2014 in visitor 
mode, which allowed us to use the blue-sensitive E2V CCD. We 
chose the 1200B grism covering the -3600-5200 A range and a 
variable slit width of 1" for June 4 and 0.7" for June 5, yielding 
a nominal wavelength resolution of AT — 3.1 A and 2.3 A for 
the two nights respectively. The slit width was adjusted to limit 
slit losses due to the poor seeing ( 1.5-2") on the first night and 
kept at the originally desired value on the second night (seeing 
of 0.8-1.2"). We employed a shutter time of 600 s per exposure 
and set up the observing sequence so as to alternate between two 
different retarder waveplate settings of -45° and +45°. On the 
first night we obtained 16 exposures (8 at each of the retarder 
plate settings), while on the second night we managed to collect 
30 spectra (15 pairs at perpendicular retarder plate settings). In 
total, we gathered nearly 8 hours of useful data on-target, the aim 
being to reach a S/N>2000/A for this relatively faint (B -13.0) 
target. We also observed the bright spectropolarimetric standard 
HD 94660 using the exact same instrument setup as a sanity 
check. 

Data reduction was done closely following the procedure de¬ 
veloped and described by [Bagnulo et al. ( 2012 ;). Simply put, we 
calculated two profiles: the circular polarisation normalised to 
the intensity Py — V/I (where V is the Stokes parameter measur¬ 
ing the circular polarization and I is the usual unpolarized inten¬ 
sity), and the null profile Ny, which is representative of the noise 
of Py. We obtained the ’rectified’ version of these quantities by 
combining the series of exposure pairs ta ken at perpendic ular re - 
tarder waveplate settings using Eq. (3) of lBagnulo et al.l (120121) . 
The results are illustrated in Fig.|T|for the second night of obser¬ 
vations (the result for the first night is qualitatively the same, but 
the S/N is lower). The top panel shows the total intensity (i.e. the 
time-averaged spectrum) together with Py and Ny as a function 
of wavelength. Even more interesting are the diagnostic plots in 
the lower panels, as is explained in what follows. The maximum 
S/N in the spectrum is -2500 per A. 

Under the so called “weak-field” approximation, the inten¬ 
sity of circular polarisation (due to the Zeeman effect) across a 
spectral line can be expressed as: 

Pv = -ge«C z A 2 ~( B z > (1) 

where g e ff is the effective Lande factor, A is the wavelength, C- 
is a constant and (B z ) is the mean longitudinal magnetic field 
expressed in G. The longitudinal field is the integrated line-of- 
sight component of the magnetic field. As can be deduced from 
Eq. (1), the latter is obtained as the slope of a linear regression 
of Vy/ versus the quantity -g e ff C- A 2 (1 //) d//d A. which is pre¬ 
cisely what is plotted in the bottom left-hand panel of Fig JT) 
Note that for the hydrogen lines of a star with 7j.fr -35,000 K 
and a dipolar field, the “weak-field” approximation is formally 
valid for field strengths up to 10-20 kG (see, e.g. Eq. 1 from 
iBagnulo et al.lll995h . 

A visual inspection of P v yields no evidence for the presence 
of a strong magnetic field. Using the least-squares technique dis¬ 
cussed in [Bagnulo et al.l (l2002h and applying it to all spectral 
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Fig. 1. Top: The total intensity FORS2 spectrum I (black) together with the normalised circular polarisation Py (red) and the Null 
profile Ny (blue) obtained by combining all data taken on 2014-06-05 for LS IV-14°116 . The null profile is offset by -0.5% for 
display purposes, the lcr errors expected from photon noise being represented by the light blue error bars. The green (magenta) 
marked wavelength ranges indicate the portions of the spectrum used for the analysis based on the H (metal) lines respectively. 
Bottom: Py (left) and Ny plotted as a function of the quantity indicated on the x-axis. The slope of Py in this diagram is proportional 
to the mean longitudinal magnetic field (B z ). The plots refer to the results obtained from using the H lines (green portions of the 
wavelength range in the top plot) and the metal lines (magenta portions) as indicated. 


lines we measure ( B z ) = 9+93 G, a null detection within the er¬ 
ror bar of about 100 G. From the H lines alone we derive ( B z ) 
= -180+205 G and from the He and metal lines we find (/?,) = 
97+105. The observations obtained during the first night also re¬ 
sulted in a null detection, although with a larger uncertainty due 
to the smaller S/N ratio of the spectra. As expected, null detec¬ 
tions were also obtained from the null profiles, confirming the 
reliability of our results. However, in Fig. 1 both the P v and N v 
profiles show more scatter around zero than would be expected 
from photon noise alone (i.e. the scatter in the profiles is larger 
than the light blue 1-cr error bars). A deeper inspection of our 
data reveals that spectra obtained during an observing series are 
slightly offset from each other. If interpreted in terms of radial 
velocity, the observed shifts would suggest that during the ob¬ 
serving series, the star has changed its radial velocity by about 40 
km/s, which is probably unrealistic and most likely due to instru¬ 
mental effects. We note that during the observing series the in¬ 
strument rotated by about 120° while the airmass changed from 


about 1.5 to 1.1. In these conditions, tests performed in imaging 
mode show that a differential shift of about 1 pixel is to be ex¬ 
pected (see Fig. 2.4 of the FORS user manual); the presence of 
a grism and polarimetric optics will presumably strengthen this 
effect. We conclude that the observed shift is probably due to to 
small instrument flexures, which in turns result in Stokes pro¬ 
files noisier than expected from photon noise. For more detailed 
information see lBagnulo et all (120131) . 

From the analysis above we claim a null detection of a mag¬ 
netic field down to a 3 cr threshold of 300 G. Note that, formally, 
this applies only to the mean longitudinal magnetic field (i.e. the 
magnetic field averaged over the visible disk in the direction of 
the line of sight) at the time of observation. The possibility re¬ 
mains that LS IV-14°116 in fact does have a strong dipolar field 
but that its longitudinal component was not detectable because 
we observed the magnetic field equator-on. Since a detailed anal¬ 
ysis of the pulsation periods of LS IV-14° 116 dGreen et alJ2011h 
showed no indication of rotational splitting, we can assume that 
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Fig. 2. Left: Fit to the Balmer and He lines in the FORS2 spectrum of LS IV- 14° 116 . The thin black line is the observed spectrum 
and the thick red line represents the best-fit model atmosphere. The quoted uncertainties refer only to the formal fitting errors, and 
certainly underestimate the true errors. Right: The equivalent fit to the UVES spectrum. 


our target is a slow rotator with a rotation period longer than a 
month or so. Therefore, our spectropolarimetric measurements 
over two consecutive nights would have been taken at a very 
similar rotational phase. Assuming a limb darkening coefficient 
of 1, the mean longitudinal field of a dipole is related to the field 


strength at the magnetic pole B p by 

(B z ) = 0.4B p cos £ (2) 

where i is the angle of the dipol ar axis with respect to the line of 
sight (see Eqs. 2 and 3 of Bagnulo et al.li'j 995 ). If we observed 
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the magnetic field pole-on we would then (at 3 cr) have been sen¬ 
sitive to polar magnetic field strengths down to ~ 750 G, while 
at 7=45° we would have barely detected a field with B p ~ 1 kG. 
A dipole magnetic field with B p > 5 kG would then have had 
a 85% chance, one with B p > 7.5 kG a 90% chance of being 
detected from our measurements. 

Another issue is that the technique used is sensitive mainly to 
the longitudinal co mpone nt of a dipolar field. According to Eq. 
5 of lLandolfi et al.l d 19981) . and again assuming a limb-darkening 
coefficient = 1, the contribution to the longitudinal field of any 
quadrupolar component is at best about six times smaller than 
the contribution due to the dipolar component. This implies that 
a quadrupolar field would have to be very strong indeed (at the 
50 kG level or higher) t o be detected by our me asurements. 
According to the work of lMathvs & Hubrigl (120061) there is an 
empirical correlation between a dipole and a quadrupole mag¬ 
netic field (which can be measured on the basis of line profiles) 
for roAp stars, i.e. stars with larger quadratic fields also have 
larger dipole fields. However, it is completely unclear whether 
such a correlation would hold for hot subdwarfs. The authors 
also point out that for a magnetic field to be detectable the organ¬ 
isation has to be large-scale enough for contributions of various 
parts of the stellar surface to not cancel out, and the geometry of 
the magnetic field along the line-of-sight has to be favourable. 

The strong magnetic fields so far measured in white dwarfs 
tend to be dipole fields with strengths on the order of several 
MG, with weaker fields at the tens of kG level also having 
been detected in a small numb er of these stars dLandstreet et alJ 
1201 2bt lKawka& Vennesll2012l). Similarly, the recently discov¬ 
ered magnetic He-sdO star dHeber et al.! 1201 31) has a magnetic 
fields strong enough to induce clear Zeeman splitting (hundreds 
of kG). We would have been extremely unlucky to have missed 
anything remotely comparable to these in LS IV-14°116 . 

Our measurements for LS IV-14°116 can be contrasted with 
those we obtained for HD 94600 , a b right Ap star with a well- 
measured magnetic field (e.g.. lLandstreet et al.ll2014l) . Here, the 
exact same least-squares fitting technique yields an integrated 

line-of-sight magnetic field of (B z ) -2200 G. This is in line 

with predictions for the time of our observations (MJD~56813) 
assuming the candidate rotational period of the star of 2800 d 
(see Lig. 5 of Landstreet et al. 2014). Thus, the instrument setup 
and technique yielding the null detection reported above for LS 
IV-14°116 appear sound. 



45000. 40000. 35000. 30000. 25000. 20000. 
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Fig. 3. We show the location of the Helium Main Sequence (solid 
line) in the log g - 74ir diagram for different stellar masses as 
indicated. The Zero-Age EHB (right-hand dotted line) and the 
Terminal-Age EHB (left-hand dotted line) are shown for differ¬ 
ent hydrogen envelope masses based on the Montreal 3rd gen¬ 
eration evolutionary models. Superimposed on this is the loca¬ 
tion of LS IV-14°116 derived from the LORS2 (left) and UVES 
(right) spectra respectively. Both estimates place LS IV-14°116 
firmly on the EHB; a star that has not yet settled on the EHB 
would be found at similar temperatures_but significantly lower 
surface gravities (cf. Fig 1. [Miller Bertolami et al.ll201 ll) . 


spectra that were originally obtained with the aim of detecting 
pulsational radial velocity variations, and the details of this are 
published elsewhere (Jeffery et al. 2014, MNRAS, submitted). 
For the purposes of the present study, we focus exclusively on 
the time-averaged absolute radial velocities and the combined 
spectrum. 


3. Spectroscopic analysis 

While the main goal of the FORS2 observations proposed for 
LS IV-14° 116 was the measurement of a possible magnetic field, 
the total intensity (Stokes 7) spectra are also of interest. In par¬ 
ticular, we planned to use the combined time-averaged spec¬ 
trum obtained to determine the atmospheric parameters, and 
hoped to exploit the time-series of individual spectra for ra¬ 
dial velocity (RV) measurements tracing the stellar pulsations. 
The latter study eventually proved fruitless due to the ~40 km/s 
RV drift over the course of the nightly observations discussed 
above. Attempts to correct for the instrumental drift yielded 
RV curves too noisy to detect pulsations at the expected few 
km/s level. Quite interestingly, we did note a strong absolute 
RV shift of around -150 km/s, which was either overlooked or 
dismissed as unimportant in previous studies of LS IV-14°11 6 
(lAhmad & Jeffervll2005l : IGreen et~ahll201 it iNaslim et al.ll201 lb . 
We were able to confirm this blueshift from a re-analysis of high 
resolution UVES spectra downloaded from the ESO archive. 
The UVES dataset comprises a 4-hour time-series of around 30 


3.1. Atmospheric parameters 

The atmospheric parameters so far derived for LS IV-14° 116 
by tw o independent groups (INaslim et alJ 1201 lb IGreen etafl 
1201 ll ) are somewhat at odds with each other, in particular as 
far as the surface gravity derived is concerned. On the basis of 
their high resolution echelle spectra, Naslim et ak| ( 201 ll) mea- 
sured T e s =34,500+500 K and log g=5. 6+0.2, while Green et akl 
(1201 lit report T eff =34,950+250 K and log g=5.93+0.04. We 
made use of the combined spectrum from both our FORS to¬ 
tal intensity spectra and the UVES archival data to re-evaluate 
the atmospheric parameters of LS IV-14° 116 . 

The FORS2 combined spectrum was computed using our in¬ 
dividual spectra obtained on June 5 only, i.e. those with the bet¬ 
ter spectral resolution of A/l ~ 2.3 A. This yielded a high qual¬ 
ity time-averaged spectrum covering the 3650-5100 A range, in¬ 
cluding the Balmer lines from H/3 upward as well as 8 strong 
Helium lines. These were fit using our standard grids of NLTE 
model atmospheres and synthetic spectra (see e.g. [Brassard et al.l 
1201 Ol for details) incorporating a fixed metallicity inspired by 
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Fig. 4. U-V diagram showing the position of LS IV-14°116 (thick error bars) as compared to 3-cr thin disk and 3-cr thick disk 
contours. Stars located outside these contours have kinematic properties of the halo population. The points inside the contours 
correspond to the reference white dwarf sample of lPauli et al.l ( 2006 ). These authors performed a detailed kinematic analysis of 398 
DA white dwarfs using the same method as applied here. 


the abundance analysis of iBlanchette et al~ l200 8t). Note that 
these are the same models as used bv lGreen et alJ (l20lfl) in their 
analysis, and that the metal abundances used are representative 
of hot subdwarfs in general, but were not tailored to the ex¬ 
tremely peculiar abundance profile of LS IV-14°116 . The com¬ 
bined UVES spectrum was analysed using the same models, but 
this time convolved to the appropriate 0.1 A wavelength reso¬ 
lution. Here, we combined the blue part of the pipeline-reduced 
individual spectra as available from the ESO Phase 3 data prod¬ 
ucts archive. 

The best-fit synthetic spectra are shown together with the ob¬ 
served FORS2 and UVES spectra and the inferred atmospheric 
parameters in Fig. [2] While the values inferred are not compati¬ 
ble within the formal errors (which are well known to underesti¬ 
mate the true uncertainties), they both yield a value of log g close 
to that inferred bv lGreen et al.l ( 201 ll ). Given the errors almost 
certainly introduced when combining the echelle orders and fix¬ 
ing the continuum of the UVES spectra, we believe the FORS2 
values to be more reliable. Therefore, we adopt 7Lrr=35,150 
K and log #=5.88 in what follows. These values place LS IV- 
14° 116 firmly on the He-burning main sequence, as can be seen 
in Fig. 0 


3.2. Kinematics 

We searched the series of archival UVES spectra for RV vari¬ 
ations by measuring the Doppler shift of the spectral lines 
with respect to their rest wavelength. This was achieved by fit- 
ting a set of math ematical functions with the FITSB2 routine 
(tNaniwotzki et al.ll2004 ) to the He I lines at 4026.2, 4120.9 and 
4387.9 A. We used a polynomial to match the continuum, while 
the line wings were fit with a Lorentz profile and the line core 
with a Gaussian profile. Heliocentric corrections were then ap¬ 
plied to the times and velocities. 

Averaging the RV measurements from the individual UVES 
spectra we infer an absolute average radial velocity shift of 


v r = -149.1+2.1 km/s. This is in line with the ~ -150 km/s 
shift seen in our FORS2 spectra. Note also that the UVES spectra 
do not show a radial velocity drift comparable to that found for 
the FORS data, confirming the instrumental nature of that effect. 
Upon reinspection, the lower spectral resolution (AT = 9 A) Bok 
telescope spectra used for the spectral analysis of iGreen et al.l 
d201 ll) show a compatible wavelength shift. There can then be 
no doubt that the blueshift detected is real. Since single subd¬ 
warf B stars in the Galactic field tend to have significantly lower 
radial velocities (|RV| < 40 km/s, Betsy Green, 2014, private 
communication) we found this highly intriguing and decided to 
run a kinematic analysis. 

The velocity components of LS IV-14°116 were com¬ 
puted based on its coordinates (a 2000 = 20h57m38.8s, (Looo = 
-14°25'47" from Simbad), radial velocity (determined from the 
UVES spectra above), proper motion and distance. All the pa¬ 
rameters were varied within the errors adopted. The distance 
was estimated to be d_= 0.44+0.04 kpc following the method 
described by iRamspeck et al.l d2001l) and assuming the values 
of log g and 7/n- as inferred from our FORS2 spectra shown in 
the left panel of Fig. [2] (but increasing the errors to more real¬ 
istic values of 0.05 and 500 K respectively), a canonical mass 
of 0.47 M q , a magnitude of V =_ 13.03 from Simbad, and us¬ 
ing the extinction maps of ISchlaflv & Finkbeineil (f2011 ). The 
proper motion was taken from the PPMXL catalogue via Vizier 
to be /To-cosa = 8.2+1.8 mas/yr and fig = -130.6+1.8 mas/yr. 
Since these observed parameters are in the heliocentric refer¬ 
ence frame, they need to be supplemented by Solar velocity and 
position information in order to be transferred to the Galactic 
rest frame. Here we assumed the distance of the Sun from the 
Galactic centre to be 8.4 kpc, its motion relative to the Local 
Standard of Rest (LSR) to have component s 1^0= 11.1 km/s, 
v v O= 12.24 km/s, v z O=7.25 km/s dSchonrich et all201(j) . and the 
velocity of the LSR to be Vi sr = 242 km/s, as predicted by Model 
I of llrrgang et al.l d2013i) . 
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The position and kinematic properties of LS IV-14° 116 were 
then computed following the equations given in Appendix A. 
We derived Galactic coordinates x = -8.10+0.03 kpc, y = 
0.20+0.02 kpc, z = -0.25+0.02 kpc and a Galactic restframe 
velocity of 60.39+19.37 km/s, with cartesian velocity compo¬ 
nents v x = -14.66+7.65 km/s, v y = -55.23+22.28 km/s and 
v z =-10.21+9.84 km/s. This translates to a Galactic radial veloc¬ 
ity component 17=13.23+8.28 km/s and a Galactic rotational ve¬ 
locity component V = -55.56+22.13 km/s. The negative value 
of V immediately indicates a retrograde orbit, i.e. LS IV-14°116 
is orbiting the Galactic centre in the opposite direction to the 
disk, implying it to be a halo star. Fig. |4] illustrates this quite 
nicely: while the stars of the reference sample fall mostly within 
the 3cr contours expected for the thin or thick disk population 
in the U - V plane, LS IV-14°116 lies well outside these con¬ 
tours. It must therefore be a halo star that by chance is relatively 
close by. For completeness, it should be noted that using the at¬ 
mospheric parameters and in particular the lower surface gravity 
of iNaslim et al . (1201 ll) for the distance estimate gives an even 
lower value of V ~ -150 km/s, i.e. the qualitative conclusion 
that LS IV-14°116 is a halo star does not change. 

4. Conclusion 

The conclusions we can draw from the study presented here are: 

- LS IV-14° 116 does not have a strong magnetic field. 
Formally, the FORS2 spectro-polarimetric measurements 
exclude a mean longitudinal magnetic field higher than 300 
G. 

- The atmospheric parameters derived from the combined 
FORS2 total intensity spectrum are consistent with LS IV- 
14° 116 residing on the Helium main sequence. 

- Based on a kinematic analysis, LS IV-14°116 appears to be¬ 
long to the halo population of our Galaxy. 

The fact that LS IV-14° 116 does not have a discernable mag¬ 
netic field means that we need to find alternative scenarios to 
explain the pulsations as well as the extremely peculiar abun¬ 
dance pattern. One currently popular idea is that LS IV-14°116 
is a pre-Helium Main Sequence (pre-He-MS) star that is still 
in the process of settling onto the Extreme Horizontal Branch. 
In this scenario, the observed abundance pattern can be qualita¬ 
tively understood in terms of an atmosphere not yet sorted by 
diffusion. Diffusion, if unimpeded by other processes, will drain 
the atmosphere of a hot subdwarf of metals (and even Helium if 
any Hydrogen is present) on a time-scale of <10 4 years, negligi¬ 
ble compared to the 10 8 -year time-scale of the He-core-burning 
phase (see, e.g. Fig. 11 of lLatour et al]l20l4 ). It is clear that in 
order for He-rich sdBs to exist in the numbers observed, diffu¬ 
sion must be slowed down, possibly by internal turbulence or a 
weak stellar wind, however the details are far from being under¬ 
stood. 

The pre-He-MS scenario is supported by the qualitative 
modelling of the pulsations observed in LS IV-14°116 in terms 
of an e-mechanism acting in the He-burning shells that ap¬ 
pear before the_star settles in the core He-burning phase 
dMiller Bertolami et al.l 12014 ). Depending on the exact model 
parameters used, periods overlapping those observed in LS IV- 
14° 116 can be excited for models in the right temperature 
range around 35,000 K. However, maximum instability occurs at 
logg ~5.2-5.7 and pulsations cease for models with logg > 5.8 
since the e-mechanism is active only during the He-burning sub¬ 
flashes following the primary He-core flash, which stop once the 


star settles on the He-MS. According to the atmospheric param¬ 
eters derived above, LS IV-14°116 resides on the He-MS with 
log g ~5.9. 

There are then two options: either LS IV-14°116 has already 
settled on the He-MS, or our atmospheric parameters are wrong. 
While we have convincing reasons to believe our atmospheric 
parameters to be rather accurate for H-sdB stars in the 29,000- 
36,000 K range thanks to frequent cross-checks with asteroseis- 
mic modelling, it is possible that the atmospheric parameters 
inferred for LS IV-14°116 will change slightly when including 
metals appropriate for the peculiar abundance profile observed 
in the model atmospheres. Note however that the assumed metal 
abundance affects mostly the derived value of T e g and has lit¬ 
tle impact on log g. A more precise abundance analysis of LS 
IV-14°116 is currently planned by another group on the basis of 
UV-spectroscopy scheduled with HST, and will hopefully pro¬ 
vide new insights on this. 

If LS IV-14°116 is confirmed as a bona fide helium core¬ 
burning star, the driving of the pulsations and the unusual sur¬ 
face composition remain an enigma. We believe that the halo 
membership may be of importance in this context. Indeed, the at¬ 
mospheric properties of halo and globular cluster hot subdwarfs 
show systematic differences compared to those in the Galactic 
disk in terms of their relative distributio n in logg-T e ff space (e.g. 
lLatour et al.ll2014t iNemeth et alJl2012h . Moreover, there is first 
observational evidence that intermediate He-rich stars in general 
belong to the halo (R Nemeth, 2014, private communication). 
The pulsational properties of hot subdwarfs in some globular 
clusters and the field also appear to be different. For instance, 
the rapid sdB pulsators of the field population have not yet been 
discovered in a_globular cluster despite systematic searches (e.g. 
iRandall et alJl201 ll) . and conversely, the fast sdO pulsators re¬ 
cently found in co Cen do no t app ear to have counterparts among 
field sdO stars dJohnson et al .120141) . It is not at all clear why this 
is the case, but must reflect differences in the internal structure 
of the stars in question. 

Presumably, the metallicity of the subdwarf progenitor popu¬ 
lation plays an important role. Apart from being metal-poor, the 
red giant progenitors of sdB stars in the halo are old compared 
to those in the Galactic disk and have a relatively low mass of 
around 0.8 M 0 . The progenitor mass is in turn predicted to influ¬ 
ence the core mass of the sdB dHan et al.l2002i) . It was also sug¬ 
gested that different formation channels for sdB stars dominate 
depending on the age of the population, the white dwarf m erger 
scenario becoming more important for older populations (iHanl 
120081') , Besides updated pulsation models including variable en¬ 
velope compositions, a better observational characterisation of 
the different hot subdwarf populations is needed to address these 
issues. It would be particularly interesting to run kinematic anal¬ 
yses for all well-studied subdwarfs with unusually large absolute 
radial velocities and/or proper motions to assess whether they are 
in fact halo or Galactic disk stars. Statistically significant sam¬ 
ples for these populations could then be subjected to detailed at¬ 
mospheric analyses and contrasted with each other. Ultimately, 
the asteroseismic exploitation of halo sdB pulsators may be used 
to compare the masses and structural compositions of these stars 
with those of their Galactic disk counterparts. 
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Appendix A: Kinematic calculations 

We calculated the kinematic properties of LS IV-14°116 based 
on the input parameters given in section 3.2 using the equations 
detailed below. To begin with, the polar coordinates of the target 
are converted to cartesian coordinates using 


C = 


d ■ 


cos a cos 6' 
sin a cos 6 
sin 5 , 


where d is the distance of the star (from the Sun) and a and 6 are 
the right ascension and declination respectively. This cartesian 
coordinate system is then rotated and shifted to x,y, z coordi¬ 
nates in the Galactic coordinate system, defined by the Galactic 
centre being at (0,0,0), the Sun being at x = -8.4 kpc (i.e. its as¬ 
sumed distance from the Galactic centre), and the positive z-axis 
pointing towards the Galactic North pole, using 


'x' 



' ~dQ_ cc ' 

y 

= M 

C + 

0 

, z, 



l 0 J 


where c/©_cc is the distance from the Sun to the Galactic centre 
and 


M = 


V T 

vf 

V T 

V NGP 


\ 
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Vgc 


' COS »GC COS 6 gc ' 
sin or GC cos d GC 
, sin S G c 


Vngp = 


' COS O-NC.R COS <5n G P 

sin qtn G p cos <5n G p 
, sin 5 ngp i 


V y - -(V G c x Vngp) 

are the rotation matrix transformation s app lied to the x, z, y axes 
respectively (iJohnson & Soderblomli 1987i) . We performed these 
transformations assuming Galactic Centre coordinates a GG = 
17:45:37.224, <5 GC = -28:56:10.23 and Galactic North Pole co¬ 
ordinates a£NP_^_J^t5J_i26.282, <5 G np = 27:07:42.01 as deter¬ 
mined by Rei d & Brunthaletl (120041) . 

For the computation of the velocities, the individual ob¬ 
served velocity components are converted to cartesian coordi¬ 
nates separately and subsequently added to give the total veloc¬ 
ity vc in the cartesian system: 


v rc — v r 


' cos a cos 5 ' 
sin a cos 5 
, sin 6 


v s = ns ■ d ■ 


- cos a sin 6 ' 

- sin a sin 6 

cos 1 5 


v a 


ju a cos d ■ d ■ 


- sin a' 
cos a 


0 ) 
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\’c = v rc + Vs + V a 

where v r is the measured radial velocity, // () is the observed 
proper motion in declination, and cos 6 is the observed proper 
motion in right ascension. The resulting cartesian velocity com¬ 
ponents are then transformed into velocity components v x , v y , v z 
in the Galactic coordinate system with 


' V, 'I 



Go ' 

Lv 

= M • 

■ v c + 

Vy© + VUr 

v z J 



K VjQ / 


where v/ jr is the velocity of the Local Standard of Rest and 
I',©, VyO, v z q are the velocity components of the Sun in the 
Galactic reference frame. The total Galactic restframe velocity 
is then given by 


Vgrf — yv x 2 + V-y 2 + V z 2 

The Galactic radial velocity and rotational velocity components 
of the star are then respectively computed as 


U 

V 


xv x + yv y 
-y/x 2 + y 2 


xvy - yv x 

yjx 2 + y 2 
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